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Abstract 

A (linear) nonsingular solution for the edge dislocation in the translational 
gauge theory of defects is presented. The stress function method is used and 
a modified stress function is obtained. All field quantities are globally defined 
and the solution agrees with the classical solution for the edge dislocation in 
the far field. The components of the stress, strain, distortion and displacement 
field are also defined in the dislocation core region and they have no singularity 
there. The dislocation density, moment and couple stress for an edge dislocation 
are calculated. The solutions for the stress and strain field obtained here are in 
agreement with those obtained by Gutkin and Aifantis through an analysis of 
the edge dislocation in the strain gradient elasticity. Additionally, the relation 
between the gauge theory and Eringen's so-called nonlocal theory of dislocations 
is given. 
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1 Introduction 

The gauge theory of dislocations [1-22] is a very attractive field of research. The 
theory of dislocations is in these theories considered as a translational gauge theory 
(T(3)-gauge theory). An important goal is to find modified solutions of the stress 
and strain field of screw and edge dislocations in this framework because the "classi- 
cal" elasticity of dislocations contains singularities at the dislocation line. Therefore, 
the classical solutions cannot be applied within the dislocation core region. This 
is unfortunate since the core region is the most important from the point of dislo- 
cation and failure mechanics. The difference in these gauge models is the form of 
the T(3)-gauge Lagrangian. Edelen et al. [6-8, 10] used the simplest form quadratic 
in the T(3)-gauge field strength. In [1-4, 13] the gauge Lagrangian was proposed 
in the three-dimensional Hilbert-Einstein form. The most general dislocation gauge 
Lagrangian for an isotropic and anisotropic material was proposed in [1,2]. 



A correct gauge theoretical solution of the stress field of a screw dislocation is 
given in [2-4,9, 10]. This solution has no singularity at the dislocation line and the 
far field stress is that of a classical screw dislocation. 

Up to now, the situation is less satisfactory for the edge dislocation. In the Edelen 
model [6-8] the far field of the stress <j zz is different from the standard one. In 
fact, the far field stress a zz is equal to v~ 1 {a xx + a yy ) instead of v(a xx + a yy ) (see 
also [4]). Unfortunately, no explicit gauge solution for the edge dislocation was given 
in [6-8]. Moreover, the structure of the gauge Lagrangian used by Edelen et al. 
requires, in principle, an asymmetric force stress [2]. Malyshev [4] used the Hilbert- 
Einstein Lagrangian as a gauge Lagrangian and found the stress field of a straight edge 
dislocation with a modified asymptotic behaviour which differs from the classical edge 
dislocation. In fact, the stress field has no singularity at the dislocation line and has 
an oscillatory behaviour with slowly decreasing amplitude far from the dislocation 
core. His solution can be considered as a complicated gauge defect with modified 
asymptotics instead of an elementary defect [4]. 

On the other hand, there are other non-standard continuum models of defects, 
e.g., the Peierls-Nabarro model [23-26], the nonlocal continuum model [27-31] and 
the strain gradient elasticity [32-38]. The force stresses calculated in these theories 
are symmetric even in the core. In the case of a screw dislocation, the solution of the 
force stresses in nonlocal elasticity [28, 30, 31] and in gradient elasticity [37, 38] agrees 
with the gauge theoretical one [2-4,9, 10]. Additionally, the strain fields calculated 
in gradient theory [35,37,38] and in gauge theory [2,3,9,10] coincide. In the case 
of an edge dislocation, the solution in nonlocal elasticity was given in an integral 
form by means of a special two-dimensional nonlocal kernel [30]. Using another two- 
dimensional nonlocal kernel Eringcn [27] gave a different solution of the straight edge 
dislocation. The dislocation stress field in gradient elasticity [36-38] was obtained 
in a closed analytical form by using the Fourier transform method, it is equal to 
zero at the dislocation line and beyond the dislocation core the classical and the 
gradient solutions coincide. In a quasi-continuum model [31, 39-42] a modified stress 
field of an edge and a screw dislocation was obtained, also. In this approach the 
stress field has no divergence at the dislocation line and at large distances it contains 
a decreasing oscillatory contribution. Dislocations have also been considered in a 
Cosserat media (see, e.g., [43-45]). The force stresses of dislocations calculated within 
the Cosserat theory are, in general, asymmetric. The solutions of edge and screw 
dislocations [46-52] found for the elastic strain and stress fields differ from the classical 
solutions but still posses singularities at the dislocation line. 

A common feature of gauge theory, nonlocal elasticity, strain gradient theory and 
theory of Cosserat media of defects is that a characteristic length scale enters the 
constitutive laws. In the gauge theory of defects and in the theory of Cosserat and 
multipolar media this length is a material property which carries with it all of the 
difference between solutions with or without moment stresses. Its influence might be 
become important as dimensions of the body or wavelengths diminish to the order 
of the characteristic inner length. By means of this characteristic length scale it is 
possible to define the dislocation core radius in a straightforward manner [2] . 

The question arises: Can we find a gauge theoretical solution for the edge dislo- 
cation which has the correct far field stress and is related to any solution of another 
framework? The aim of this paper is to find an answer to this question, that means 
to construct a correct gauge theoretical solution of a straight edge dislocation. Addi- 
tionally, the relation to the nonlocal theory and strain gradient elasticity of an edge 
dislocation will be given. Here we restrict our considerations to the case of linear 
dislocation theory. 
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2 The translational gauge invariant formulation of 
dislocation theory 

For the undeformed reference state of the body the indices i = 1, 2, 3 are used. The 
final coordinates of the deformed state are labelled by the indices a = 1, 2, 3. In pres- 
ence of defects the deformed state of the material is characterized by nonholonomic 
coordinates. In linear theory we may use S a l , Sij and 5 a b to raise and lower the indices. 

In elasticity theory the state of an elastic body is not changed under a rigid 
translation and rotation (global gauge transformations) . In elasticity the (compatible) 
distortion is given as a proper gradient 

P% = d t u a . (2.1) 

The distortion tensor is dimensionless. The basic idea of gauge theory of defects is now 
that the distortion has to be invariant not only under the global gauge transformation 
but also under if the gauge group is applied locally. In this way, one obtains the 
incompatible distortion in the framework of translational gauge theory of dislocations 
as follows 

P*. = d t u a + (2.2) 

It is locally translational invariant (gauge invariant) due to 

u a ^u a +T a (x), 0«. — 0«. - ^(s), (2.3) 

where r°(x) are local translations [2]. Here u a is the displacement field from the 
undeformed to the deformed configuration and (\> a i is the proper incompatible part of 
the distortion. The incompatible distortion (2.2) can be understood as the (minimal) 
replacement of the compatible distortion (2.1) in the T(3)-gauge theory 

d t u a — > 8 t u a + a 4 . (2.4) 

Then, the (incompatible) strain tensor turns out to be 1 

Eij = ^ai/^+^i) =%). (2.5) 

The skew-symmetric part of the distortion defines the elastic rotation tensor 

w« = ^(*«i/9°j - = %]> ( 2 - 6 ) 
and the elastic rotation vector 

i^i^—e^Saj^. (2.7) 
The elastic strain energy of an isotropic media is defined by 

W=^n (s ik S jl + 6 il 5 jk + 6 ij 5 kl ^j EijE kU (2.8) 

where fi is the shear modulus and v is Poisson's ratio. It is now invariant under local 
translations. Obviously, no elastic rotation enters the relation (2.8) and, therefore, 

1 We will be using the notations = ^(Aij + Aji) and A^ = ^(Aij — Aji). 



3 



the elastic rotation has no elastic response quantity. The force stress can be derived 
from the strain energy 

,«=™1 = 2m(^ + ^^). (2.9) 

It is the elastic response quantity to the strain tensor. Formula (2.9) expresses the 
(incompatible) stress-strain relation (generalized Hooke's law). 

It is well known that nontrivial traction boundary problems in the variational 
formulation of the gauge theory of defects can be formulated by means of a so-called 
null Lagrangian [8, 10,53]. When the null Lagrangian is added to the Lagrangian of 
elasticity, it does not change the Eulcr-Lagrange equation (force equilibrium) because 
the associated Euler-Lagrange equation, dja J = 0, is identically satisfied. After 
minimal replacement, diU a — > /3", the null Lagrangian 

W hg = dj{aju a ) = {djaj)u a + aJd jU a — aj^, (2.10) 

give rise to a background stress tensor a J which can be considered as the nucleation 
field in the gauge theory of defects [10]. 

As usual we use the translational field strength (torsion) to define the dislocation 
density tensor [1-3, 16, 17, 54-57] 

T%=d i p a j -d j f3 a i , (2.11) 

which has the dimension of an inverse length. The dislocation density is a fundamental 
quantity in plasticity because the dislocation is the elementary carrier of plasticity. 
This field is a physical state quantity which, at least in principle, can be measured 
without knowing anything about the history of the body. It is possible with the help of 
high resolution transmission electron microscopes (HRTEM) to see single dislocations 
with their core configuration in crystals. Since dislocations change the energy of 
the body (crystal), they should appear in the total Lagrangian. The dislocation 
Lagrangian, which we use, contains only the translational gauge field strength and is 
given by 

C disl = ~T%H^, (2.12) 

where the moment stress H a l i is the response quantity to dislocation density. This 
means, that at all positions where the dislocation density is non zero, also localized 
moment stresses are present. It is convenient to perform an irreducible decomposition 
of the torsion, from which we can construct the most general isotropic constitutive 
law between the dislocation density and the moment stress in the following way [1,2] 

3 

flaij = $>' (J)T °«- (2-13) 

1=1 

Here a\, a-i and CI3 are new material coefficients. They have the dimension of a force. 
The irreducible decomposition of the torsion under the rotation group SO(3) reads 

T ~ Wt ■ ■ -i- ( 2 )t ■ ■ 4- ( 3 )t ■ ■ (1 14 s ! 

The tensor, the trace and the axial tensor pieces are defined by (see also [2, 58]) 

W T aij = T aij - WT aij - ^T aij (tentor), (2.15) 

{2) T aij := \ (S at T\ 3 + 5 ao T l a ) (trator), (2.16) 

{3) T alJ := 1 (T a ij + T ija + T jai ) (axitor). (2.17) 
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Then, the moment stress tensor reads 

H ™i = - 2 7«S = c i Taij + c 2 {T ija + T jai ) + c 3 (5 ai T\ + S aj T l u ), (2.18) 



with the abbreviations 

ci := ^(2ai + a 3 ), c 2 := ^ (a 3 - <n) , c 3 := ^(a 2 - <n). (2.19) 

The Euler-Lagrange equations for the Lagrangian C = £disi — W + Wb g turn out 
to be [2] 

= djaj = (force equilibrium), (2.20) 

—— = diH^ = a J (moment equilibrium), (2.21) 



with the so-called effective stress tensor 



a u a > 



(2.22) 



which drives the dislocation fields. The field equations (2.20) and (2.21) determine the 
fields /?" and T" • which represent the geometric (or kinematic) degrees of freedom. 
Note that in our approach the dislocation density is not given, a priori, as delta 
function. In this field theoretical context, we are able to calculate even the dislocation 
core. The static responses are the fields a J and . 

The L.H.S. of (2.21) can be written in terms of the distortion tensor 

d l H aiJ = Cl (Ap aj - + c 2 (d'djfca - d a d l l3 l3 + d a d*l3 Jt - A(3 ja ) (2.23) 

+ c 3 (d a d l f3 l3 - d a djFi + s aj Ap\ - d a3 d*d k f3 lk ). 

In the general case the gauge theory of dislocations for isotropic materials contains 
three new material coefficients which can be used to define three internal character- 
istic lengths. In order to obtain only one coefficient from the three ones we need a 
suitable choice for the moment stress tensor H a ij and the coefficients a\, a 2 and a 3 , 
respectively. By means of such a choice we want to find a gauge theoretical solution 
of symmetric force stresses for the edge dislocation. If we use the choice 2 

l + v ai ai ai v , nnA \ 

a2 = yyai, a 3 = — 2~. ci = y, c 2 = — -, c 3 = -ai, (2.24) 

the localized moment stress tensor reads 

Haij = y (^Taij — Tija — Tj a i + y— - (S a iT\j + 5 a jT l a )^j . (2.25) 

With the relations H a k — \^ijkH a 1 ^ and = eijka ak , and by using the conventional 
(microscopic) dislocation density tensor a a% — e l ^ k dj(i a k and the Nye [59] tensor 3 n a j 

K aj = dja — g ^°-d°^ki a aj = K ja — ^aj^ki (2.26) 



2 Note that Malyshev [4] used the so-called Einstein choice <J2 = —01 and aa = — ^y- for an 
(modified) edge dislocation. By applying the stress function method he found a gauge theoretical 
solution for the edge dislocation in terms of the Bessel and Neumann functions and the modified 
Bessel functions. His solution of an (modified) edge dislocation has no singularity at r = but differs 
asymptotically from the classical edge dislocation. In detail the symmetric stress field, the closure 
failure and the dislocation density show a decaying behaviour with oscillations. Additionally, his 
solution violates in the dislocation core region the plane-strain condition E zz = 0. 

3 Here the tensor K a j is a microscopic quantity. 
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we obtain from (2.18) the relation 



H aj = ai ( a aj - - 5 aj a\ + -^—^ a [aj] I = a x ( K ja + — — n [ja] ) . (2.27) 



a [aj] j=a t {K ja + 
Alternatively, the dislocation Lagrangian (2.12) can be written in the form 

C disl = ~a ai H ai . (2.28) 
By using (2.27) the L.H.S. of (2.21) reads 

d l B ai3 = ai j A/? (oj -) - PdjpM + ^ a Q a l %] (2.29) 
- -[^({dad, - 5 aj A)(3\ + 5 aj d t d k (3 lk - 5 a ^%)) j. 

Obviously, the strain and the rotation tensor appear in (2.29) because the distortion 
is not symmetric. It is reminded that in the so-called Hilbert-Einstcin choice only the 
strain tensor enters the moment equilibrium (see [2,3]). 
We now define the rotation gradient 4 

k%=d a p [l3] . (2.30) 

It has the dimension of reciprocal of length. The rotation gradient describes a "lattice- 
curvature". It contains two pieces. To see this we rewrite (2.30) 

hjk = {dkE l j — d 3 E ik ) + -(Ti-jk - T jki — T ki j). (2-31) 

The first one is caused by local (incompatible) clastic strain and the other one by 
dislocations (contortion) [57]. The second piece may be identified with the plastic 
"lattice-curvature" which cannot be realized without dislocations. We add the fol- 
lowing gauge invariant Lagrangian which looks like a rotation gradient Lagrangian of 
a Cosserat medium (see, e.g., [62]) 

£grad = ~\k a i^: j - (2.32) 

The response quantity to the rotation gradient is a so-called "Cosserat" couple stress 
for an isotropic material 5 (sec, e.g., [62]) 

dC 

VaM = ~ dk aM = + d 2(h[ja] + + d 3 (S ai k l M + 5 a] k l [u] ) . (2.33) 

This constitutive relation uses the three rotation gradients, which are invariant under 
the SO(3), and is analogous to the invariant constitutive relation for the torsion (2.18). 
We find 

d l Haj] = diA/? [ojl + (da - d 3 )(d*d 3 P [la] - d l d a p M ). (2.34) 

Another gauge invariant Lagrangian, which we add, is given in terms of the torsion 
and the rotation gradient 



C gIltd - disl = --dok%T^. (2.35) 



4 DeWit [60, 61] called such a rotation gradient elastic bend-twist. 

5 The relations between the constants in (2.33) and those in [62] arc: di = 2ct2, (I2 = «3, ds = — ai. 
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Equation (2.35) may be considered as a kind of an interaction Lagrangian between 
dislocations and the "lattice-curvature" . We define the response quantity 

Xaij = 8(8*13^ ) = d ° i katJ + 2 Tiaj ) ' *~ 2 ' 36 ^ 

with 

d'Xaij = d 8 l (d a (3 m + l - (dafa - djfoa) ) . (2.37) 

The coefficients do, d\, d 2 and d^ have the dimensions of forces. 

If we use the Lagrangian C = £disi + £ g rad + £ g rad-disi — W + Wbg, the moment 
equilibrium reads 

d % (H ai j + X ai j + Hi[aj\) = <Taj, (2.38) 

which can be decomposed into its symmetric and antisymmetric part. For the sym- 
metric part of the moment equilibrium we obtain the equation 

+ \aj)i) + V(aj) = 0- (2.39) 

By using the choice of the coefficient 

d ° = -jrh-y (2 - 40) 



formula (2.39) can be rewritten in terms of the strain tensor, only, according to 

'{ 



oi|A£7« - d k d {l E j)k - JL-^didj - 5 l3 A)E k k + 5 l3 d k d l E kl - d k d {i E j)k j | = 'd l3 . 

(2.41) 

This equation can be interpreted as the (static) field equation in the elastoplastic 
theory of dislocations. It is derived from the combination of several Lagrangians - 
Eqs. (2.12), (2.32), (2.35) - and the consideration of the null Lagrangian (2.10) and the 
strain energy (2.8) in combination with the special choices of material constants (2.24) 
and (2.40). The antisymmetric part of the moment equilibrium 6 

9 l (H[ aj ]i - \a\i\j] - Mi[oj]) = 0, (2.42) 

can be satisfied by the help of the choice of the coefficients 

di=0, d 2 -d 3 = y. (2.43) 

The physical interpretation of Eq. (2.42) is the following. Because we have required 
that the force stress should be symmetric, a^jj = 0, we have had to introduce the 
couple stresses A[ a |j|j] and Hi[ a j] m order to fulfill the antisymmetric moment equi- 
librium (2.42). Consequently, (2.42) expresses the moment equilibrium between the 
antisymmetric moment stresses H[ a j]i, \a\i\j] an d Mi[aj]- The fact that the moment 
stress is the specific response to the presence of dislocations and that the field equa- 
tion (2.38) expresses the moment equilibrium is a hint that the gauge theory of dis- 
locations may be connected with the theory of a constrained, anholonomic Cosserat 
media with symmetric force stress. The dislocations bring the anholonomity into the 

6 Indices which are exempt from antisymmetrisation are enclosed by vertical bars. 
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picture. A similar remark in this direction was given by Kluge [45]. It should be 
also possible to extend this theory to a more general dislocation theory in a general 
Cosserat continuum with antisymmetric force stress and 10 material constants which 
is not the aim of the present paper. However, for antisymmetric force stress and 
vanishing dislocation density the moment equilibrium equation recovers the form of 
the Cosserat equilibrium equation under zero external moments. 
Finally, we may define the hyperstress 

dC 

W = -^T) - H *' + V J + MV- (2-44) 

This hyperstress tensor is close to Maugin's so-called "Piola-Kirchhoff" hyperstress [63]. 
Such a hyperstress tensor was also used by Toupin for the theory of nonsimple elastic 
materials [64]. Thus, the dislocation theory given in this paper can be considered as 
an anholonomic higher-order gradient theory in nonsimple materials. If ai = 0, there 
are no effects of hyperstresses and couple stresses. 
By the help of the inverse of Hooke's law 

and using the equilibrium condition d J <Jij — 0, we obtain from Eq. (2.41) the inho- 
mogcneous Helmholtz equation as field equation for every component of the stress 
tensor 

(l - K- 2 A)aij = aij, k 2 = ^. (2.46) 

It is important to note that (2.46) agrees with the field equation for the stress field in 
nonlocal elasticity [28, 30, 31] and in gradient elasticity [37, 38]. The factor k _1 has the 
physical dimension of a length and therefore it defines an internal characteristic length 
(dislocation length scale). By using (2.45) and (2.46) we obtain an inhomogeneous 
Helmholtz equation for the strain fields 

(l-K- 2 A)E ij =E ij , (2.47) 

o 

where £y is the background strain tensor. Equation (2.47) is similar to the equation 
for strain of the gradient theory used by Gutkin and Aifantis [35-38], if we identify 
kT 2 with their corresponding gradient coefficient (see, e.g., equation (4) in [37]). Note 
that only for compatible strain, <j>ij = 0, one can deduce an inhomogeneous Helmholtz 
equation for the displacement from (2.47) which is analogous to the equation for the 
displacement in gradient elasticity (see equation (3) in [37]). 

The conditions on a gauge theoretical solution of dislocations in this framework 
are: (i) the stress field should have no singularity at r = 0, (ii) the far field stress 
ought to be the classical stress field. We use the condition (ii) because the classical 
stress fields are well-established and in good agreement with physical observations, 
e.g. the investigation of stresses around dislocations by means of the photoelasticity 
method (see, e.g., [65-67]). 
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3 The gauge theoretical solution of the straight edge 
dislocation 



3.1 The stress, strain and displacement field in an infinitely 
extended body 

Let us now derive the correct dislocation fields within the translational gauge theory. 
We take Cartesian coordinates, so that the z-axis is along the dislocation line and 
the z-axis is along the Burgers vector: b x — b, b v = b z = 0. The extra half plane lies 
in the plane x = 0. In order to satisfy the equilibrium condition, we use the second 
order stress function / and specialize to the plane problem of an edge dislocation by 
setting d z = 0. We are using the so-called stress function ansatz [26] 

d 2 f -d 2 f 

yyJ ^xyJ 

-dlyf d 2 x J ], (3.1) 




with 



P = vAf, (3.2) 



where A denotes the two-dimensional Laplacian d 2 x + d 2 y . The relation (3.2) comes 
from the requirement: E zz — (plane strain). Additionally, the strain is given in 
terms of the stress function as 

d 2 vv f - vAf -d%f 

E^ = ~\ -d%f dlJ-vN |. (3.3) 


We use the classical stress field of a straight edge dislocation in terms of the Airy 
stress function as the background stress 



'H - -dlx d 2 xxX . (3.4) 




The Airy stress function [68-73] 

X = -Aylnr, A = ^£— y (3.5) 

fulfills the following inhomogeneous bipotential (or biharmonic) equation 

AAx = -ATtAd y 5(r), (3.6) 

where 5(r) is the two-dimensional Dirac delta function and r 2 = x 2 +y 2 . Substituting 
(3.1) and (3.4) into (2.46) we get 

(A- k 2 )/ = K 2 Ay\nr. (3.7) 

f = -Aylnr + f (1) (3.8) 

^A - k 2 ) = 2 A d y In r. (3.9) 



We make the ansatz 



and obtain 
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Now we set 



with 



f (1) = 2Ad y f {2) (3.10) 
(A- k 2 )/ (2) =lnr (3.11) 



and we get the solution 



/ (2) = ~(lnr + ^ («r)). (3.12) 
Finally, we obtain the solution of (3.7) 

f = -A (yhir + ^ d v (inr + KoiKrj^ . (3.13) 

Then the solution of the modified stress function of an edge dislocation with Burgers 
vector b\\x reads 7 

/ = -^£^y {hr + ^(l - nrK l{ nr)) } , (3.14) 

where the first piece is the Airy stress function and K n is the modified Bessel function 
of the second kind and n = 0, 1, . . . denotes the order of this function. 

By means of Eqs. (3.1) and (3.14), the modified stress of a straight edge dislocation 
is given by 

fib y 



>xy 



{(y 2 + ^ 2 ) + ^(y 2 -^ 2 ) (3.i5) 

- 2y 2 nrK 1 (nr) - 2(y 2 - 3x 2 )K 2 (Kr)^ 

{(y 2 -- 2 )-^(y 2 -^ 2 ) (3-i6) 

- 2x 2 nrK 1 {nr) + 2(y 2 - 3x 2 )K 2 (nr)^, 

mi--)H {x2 -"' ) -^ (x '-" 2) (3 - i7) 

- 2y 2 nrK 1 {nr) + 2(x 2 - 3y 2 )K 2 ( K r)^, 

--^-^{l-nrK l{K r)}. (3.18) 
7r(l — v) r A V > 



2tt(1 - v) r 



y f/..2 „2\ 4 

"2tt(1 - v) 



jib 



ybv y 
(1-u) 

The stress a zz satisfies the condition a zz = v{o xx + cr yy ). The trace of the stress 

k 



tensor a k k produced by the edge dislocation in an isotropic medium is 



*~3t^?{ 1 -~*<"4 < 3 - 19) 

The spatial distribution of stresses of an edge dislocation near the dislocation line 



7 For an edge dislocation with the Burgers vector b\\y the modified stress function is given by: 

fib 
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-10 -5 5 10 -10 -5 5 10 

(a) (b) 




are presented in Fig. 1. 

Let us now discuss some details of the core modification of the stress holds. The 
stress fields have no artificial singularities at the core and the maximum stress occurs 
at a short distance away from the dislocation line (see Fig. 2). In fact, when r — > 0, 
we have 

^. 4 ( 3 ' 2 °) 

and thus ov, -> 0. It can be seen that the stresses have the following extreme val- 
ues: \a xx (0,y)\ ~ 0.546^2^^ at \y\ ~ 0.996k" 1 , \a yy (0,y)\ ~ 0.260*5^^ at 
|y| ~ 1.494k" 1 , \a xy (x,0)\ ~ 0.260k at |x| ~ 1.494k" 1 , and |ct zz (0,j/)| ~ 
0.399k ^-v) a ^ I S/ 1 — 1.114k" 1 . The stresses a xx , a yy and a xy are modified near 
the dislocation core (0 < r < 12k" 1 ). Note that \a xy (x, 0)| ~ 0.260/t ^n^) can 
be identified as the theoretical shear strength. The stress a zz and the trace a k are 
modified in the region: < r < 6k" 1 . Far from the dislocation line (r ^> 12k" 1 ) 
the gauge theoretical and the classical solutions of the stress of an edge dislocation 
coincide. Thus, the characteristic internal length k" 1 determines the position and 
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<Jxx(0,y) 




-0.2 
-0.4 




Figure 2: The stress components near the dislocation line: (a) a xx (0, y), (b) tr yy (0, y), 
(c) >?xy(x,0) are given in units of /x6k/[27t(1 — v)] and (d) <j zz {0,y) is given in units 
of \xbvnj\K(\ — y)]. The dashed curves represent the classical stress components. 
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the magnitude of the stress maxima. It is interesting and important to note that the 
gauge-theoretical solutions (3.15)-(3.18) agree precisely with the gradient solutions 
given by Gutkin and Aifantis [37,38] (with k~ 2 = c, c is the gradient coefficient). 
For the clastic strain of an edge dislocation we find 

E ** - - 2 -) + 2 4 + 4rM - 3- 2 ) (3.21) 



47r(l - v) r 2 



2 - ") «r/ri(«r) - ^|(y 2 - 3x 2 )^ 2 ( K r)J, 



^ /srtfifar) + A (a: 2 - 3y 2 )/f 2 (Kr) J , 



which is in agreement with the solution given by Gutkin and Aifantis [36-38] in the 
framework of strain gradient elasticity. The plane-strain condition E zz = of the 
classical dislocation theory is also valid in the dislocation core region. The compo- 
nents of the strain tensor have the following extreme values (v = 0.3): \E xx (0,y)\ ~ 
- S0SK 4^hu) at \V\ - 0.922k- 1 , \E yy {0,y)\ ~ O.OIOk^^^ at \y\ ~ 0.218k- 1 , 

at 



|£7 w (0,y)| ~ 0.054 K37( ^ at |y| ~ 4.130k- 1 , and \E xy {x,Q)\ ~ 0.260 K3 ^y 



x ~ 1.494k 1 . It is interesting to note that E yy (0, y) is much smaller than E xx (0, y) 
within the core region (see also [37,38]). The dilatation E k k reads 

It is convenient to rewrite the stress and strain fields in cylindrical coordinates. 
The stress function reads in cylindrical coordinates: 



/ = -; M& 



sin <p |r In rH — ^- ^1 — nrKi(nr)j | . (3.25) 



2tt(1 - J/) 

Here the stress function is related to the stresses by the equations 

Crr = - <9 r / + ^ d 2 ^/, CT W = <9 2 r /, (T r <p = ~d r ^- 9 v .f^ , cr Z z = K°Vr + C^v)- 

(3.26) 

In this way we find the stress of a straight edge dislocation in cylindrical coordinates 
\ib sin ip 



{i-^+^Kr)}, (3.27) 
:f l * +2*^)1, (3.28) 



2?r(l -v) r 



2tt(1 -v) r 



'ipip 



'''' S ' mlfi f l + -1^ - 2^ 2 («r) - 2«rJf 1 («r)l , (3.29) 
K z r z J 



2tt(1 -i/) r 
sin (p 
7r(l — i/) r 



— {l /••/•A; :/,•/•:}. (3.30) 
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and the trace of the stress tensor reads 

/j,b(l + v) sin ip 
7r(l — v) r 



a h = 



-jl - KrKiinr)}. (3.31) 



Eventually, the elastic strain of a straight edge dislocation is given in cylindrical 
coordinates as follows 

E rr = —7-7^ r — ((1 - 2i/) - + 2K 2 (nr) + 2vKrK 1 {nr)\ , (3.32) 



L ■•■ = -4^0 ^ " ^ + ^ " 2/V2< "' r) " ^ 1 ~" ( "' r) 



(3.34) 

The dilatation reads 

Let us now calculate the distortion of an edge dislocation. The distortion /?" is 
given in terms of the stress function (3.14): 

/ dLf-vAf -d 2 xy f + 2^ 
Pi = — I - 2 ^ 9 L/ - ^A/ | , (3.36) 



2/z 







where u> is used to express the antisymmetric part of the distortion, u> = (3[ xy ] ■ Even- 
tually, u) is determined from the conditions: 

T x xy = i- (2»d x w - (1 - v)d y Af) , (3.37) 

n y = ^ ( 2 K^ + (1 - f)»xA/) = 0. (3.38) 
We find for the elastic distortion of the straight dislocation 

47r(l - i/) r 2 r 2 K 2 r 4 v 7 



(3.39) 



y 2 



r> = Miby ? {< 3 - 2 "» - ¥ - ^ - 3 " 2 > <"°> 



- 2 ^(1 - v) + t^j nrK^nr) + ^(x 2 - iy 2 )K 2 {nr) 
" , . = -4^)?{ <1 - 2 " , + ¥ + ^^- 3 '' 2 ) < 3 ' 41 > 



2 /2 , 21 



2 ( (1 -i/) — ^2 J ktK\ («;r) - ^ (a; - 3y 2 ) J ftr 2 (^) 



2 



2 



t/ ) KrKi(nr) + ^(y 2 - 3a; 2 ) ^(k^) 
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and for the rotation 

b 



-LJ = — 



— ^{l-nrK^r)}. (3.43) 



Eq. (3.43) is in agreement with the rotation vector calculated in the linear theory 
of dislocations in the Cosserat continuum [50]. The far fields of (3.39)-(3.43) are 
identical to the classical ones given in [74]. 

We now rewrite the distortions (3.39)-(3.42) as follows 

1 / 4 2 \ 

+ TiTi r sin 2y) -=-^ - kK^kt) K 2 (nr) dr 

2(1 - v) \K z r A r J 

= du x + (jf. (3.44) 
We may interpret the field 

4> x = -—iprK (nr) dr (3.45) 

as the proper incompatible part (negative plastic distortion). This proper incom- 
patible distortion is exactly the same as (f> z of a screw dislocation (see [2,3]). In 
components it reads 

4>\ = x V K ( K r), <ff° y = -— W K {nr). (3.46) 

The appearance of such proper incompatible distortion is a typical result of the gauge 
theory of dislocations. When kT 1 — ► 0, the incompatible distortion (3.46) converts to 
<j) x x = —bx<p5(r) and (jF = —byip5(r). The compatible part of Eq. (3.44) is given in 
terms of the displacement field 



b_ 

2w \ 



|(1 - ktK 1 {kt)) V > + -^—^ sin 2^ (l - -i_ + 2K 2 {nr)^ }. (3.47) 



It is a modified displacement field (see Fig. 3a). We used tp — arctan(y/x) as a multi- 
valued field. Note that the first part in (3.47) is different from the corresponding one 
given in gradient elasticity by Gutkin and Aifantis [36,37]. One reason is that they 
used a compatible elastic distortion. It is important to bear in mind that u x is not 
a simple arctan function as in the Volterra model of a dislocation. Additionally, we 
find from the distortion 

P y = — - — 7T- A sin 2ip (1 - -5L + 2K 2 {nr)) d<p 



47r(l — v) \ \ n 2 r 2 

- ^(1 - 2v) (J- - nKi(nr)^ +cos2^ ^7^3 - kK^kt) - ^K 2 (Kr)^j ^jdr 
= du y + <j> y , (3.48) 
the incompatible part 



(3.49) 



and a single- valued displacement field 



P = -^A-^2(l - 2i/)(lnr + K (nr)) + cos2^ (l - -±j + 2K 2 (^ J 



(3.50) 
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u x (x,0)/b 
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Figure 3: Effective displacement fields [y — 0.3): u x (x,y 
The dashed curves represent the classical solution. 



+0)/b, (b) u»(a;,0)/6. 



which has no singularity (see Fig. 3b). This means that (3 V is a proper compatible 
distortion. That is the reason why the field u v agrees precisely with the correspond- 
ing formula in gradient elasticity (see [36,37]). The displacement field u v differs from 
the classical one in the region < r < 12k" 1 . In this framework the displacement 
fields (3.47) and (3.50) have no singularity. The far fields of the displacements (3.47) 
and (3.50) are identical to the classical ones given in [26,75,76]. It can be seen that 
the solution (3.47) leads to an asymmetric smoothing of the displacement profile (see 
Fig. 3a). In fact, it is asymmetric around the dislocation line. Note that the incom- 
patible distortion (3.46) also has an asymmetric form. These asymmetries would lead 
to an unexpected asymmetry of the dislocation core. From the mathematical point of 
view, the asymmetry is induced by the decomposition of the elastic distortion (2.2) 
into compatible and incompatible parts. Keep in mind that such a decomposition 
is gauge invariant (see Eq. (2.3)) and we have a local translation as an additional 
degree of freedom. This gauge invariance can be used to improve the situation. So 
far we have considered ip as a multi-valued field. On the other hand, in the defect 
theory [35-38, 74] ip is usually used as single-valued and discontinuous function. It is 
made unique by cutting the half-plane y — at x < and assuming ip to jump from 
7r to — 7r when crossing the cut. If we use the single-valued discontinuous form for p 
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u x {x,0)/b 




-10 -5 5 10 

KX 



Figure 4: Displacement u x (x,y — > +0)/b near the line of edge dislocation within 
gauge model (solid), the gradient model (small dashed) and the classical elasticity 
(dashed). We used k 2 = 1/c. 



and the local gauge transformation (2.3) with the gauge function 



b 7T 

7^ si g n (y) o Kr - K 'i( Kr )) ( 3 - 51 ) 

Z7T Z 



we obtain for the improved displacement field 



— — KrKi(nr)) + — sign(y) nrKi(nr) 



2(l-i/) r 

and for the proper incompatible distortion 
b 



»)}• 



= K 2 xK (Kr) - | sign(y)) , (3.53) 
j^yifo^r) (<p - ^ sign(y)) + 7r<5(y)(l - sign(x)[l - nrK\ (kt - )] ^ | . 



_6_ 

27 



(3.54) 



Now the displacements and the incompatible distortions are symmetric in the core 
region around the dislocation line 8 . It is interesting to compare the displacement 
field (3.52) with the displacement calculated by Gutkin and Aifantis. When y — ► 0, the 

8 Note that in the case of a screw dislocation the symmetrical displacement and the corresponding 
incompatible distortion read, respectively, 

u z = — j<^(l — KrKi(Kr)) + — sign(y)/cr\ffi(/cr)J , 



and 



~ | K 2 yK (kt) (v>—^ si S n (?/)) + n6(y)(l- sign(x)[l - KrKi(Kr)] ) J 
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Figure 5: Effective Burgers vector b x (r)/b (solid). 



Bessel function terms in (3.52) lead to the symmetric smoothing of the displacement 
profile in contrast to the abrupt jump occuring in the profile of the classical solution 
(see Fig. 4). When y — > 0, the displacement u x (x,y) calculated by Gutkin and 
Aifantis [36,38] within the gradient model has the explicit form 



The exponential term which appears in the gradient solution (3.55) leads to the 
smoothing of the displacement profile (see Fig. 4). The size of the transition zone 
is approximately 12k -1 which gives the value r c ~ 6k" 1 for the core radius. Conse- 
quently, in the gauge theory of dislocations and in gradient elasticity the dislocation 
core appears naturally. Therefore, these displacement fields may be used to model the 
dislocation core. In this way, one can compare the calculated dislocation core with 
HRTEM micrographs and related computer simulations of the core region. It would 
be very interesting to check a possible material dependence of kT 1 . 
Finally, the effective Burgers vector can be calculated as 



This effective Burgers vector b x (r) differs from the constant Burgers vector b in the 
region < r < QkT 1 (see Fig. 5). For the value of kT x = 0.25a the core radius is 
r c = 1.5a and for the value of kT 1 = 0.399a the core radius is r c — 2.4a. Note that 
the effective Burgers vector b x (r) of an edge dislocation has the same form as the 
effective Burgers vector b z (r) of a screw dislocation which is obtained in [2, 4, 9, 10]. 

3.2 The dislocation density, moment stress and core energy 

The proper incompatible part of the elastic distortion gives rise to a localized torsion 
and moment stress tensor. We find for the torsion 




(3.55) 




7 



(3.56) 




(3.57) 



(3.58) 
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which satisfies the translational Bianchi identity: dia al = 0. On can see that the 
dislocation density is short-reaching. It is interesting to note that the torsion T x xy of 
an edge dislocation has the same form as the torsion T z , of a screw dislocation which 
is given in [2-4,9,10]. In the limit as k — ► 0, the gauge theoretical result (3.58) 
converts to the classical dislocation density a xz — bS(r). Additionally, the dislocation 
density (3.58) agrees with Eringen's two-dimensional nonlocal modulus (nonlocal ker- 
nel) used in [28-31]. Consequently, the dislocation density tensor is Green's function 
of the Hclmholtz equation: 

(l- K - 2 A)a xz (r)=b5(r). (3.59) 

This torsion gives rise to a localized moment stress. The localized moment stress of 
bending type is given by the help of (2.27) as 

H xz = ,f K ( K r), H zx = J^-Ko(Kr). (3.60) 

7T(1 - V) 7T(1 - V) 

This expression is very close to the moment stress of an edge dislocation given in [77] . 
By means of the dislocation density and the moments stress we are able to calculate 
the core energy of an edge dislocation. The dislocation core energy density is given 
by 

1 nh 2 K 2 

V COIC = - a al H al = — K ( K r) 2 . (3.61) 

2 47T^(1 — V) 

Obviously, the moment stress H zx gives no contribution to the dislocation core en- 
ergy. We obtain for the dislocation core energy per unit length of the straight edge 
dislocation 

r R 

Score = / dip I r Vcoredr 



rZ7T pi 
= dip 

Jo Jo 



Lib 2 *, 2 



{K Q ( K r) 2 - K^r) 2 } 



4tt(1-i/) L J lo 



^—Al + ^^^inRf-K^KR) 2 )^, (3.62) 



4tt(1 - v) 

where R denotes the size of the solid body. In the limit R — > oo, the core energy of 
the straight edge dislocation, 

Ecmc = Mi-vy (3 ' 63) 

coincides with the core energy that is calculated for the edge dislocation in the Peierls- 
Nabarro model [26] . The core energy is different from the constant value in the region 
< kR < 3. 

The far-reaching rotation gradients read 



k* = -k* = -- 



b 



xy ' 2irr 4 



|(x 2 - y 2 )(l - KrK^nr)) - n 2 x 2 r 2 X (/cr)}, 



k y xy = ~k\ x = -^^{2(1 - nrK^r)) - K 2 r 2 K a {nr)). (3.64) 

The form of these rotation gradients is in agreement with the expressions calculated 
within the theory of Cosserat media (see [48,51]). 



19 



3.3 The straight edge dislocation in a cylinder 

Now we consider a straight edge dislocation in a cylinder of finite outer radius R with 
an outer boundary free from external forces. An advantage in our framework is that 
we do not need an inner radius because the stress and strain fields are zero at r = 
in quite natural way. If the outer boundary is to be free, it is necessary to add to the 
stress function (3.5) an additional stress function \' = ACyr 2 . Then the modified 
stress function is of the form 



/ --TZTT—J >in ; |rl:i/- ' 4^ ( 1 KrK^nr) \Cr) - CV 3 j , (3.65) 



/j,b 
2tt(1 — v) 



where the constant C is determined by the semi-classical boundary condition that the 
stresses a rr and a rip should be zero at r = R. In this way we find the stress of the 
edge dislocation in cylindrical coordinates 

*r* = ^ (l - 4-, + - ^r 2 ) , (3-66) 

Z7r(l — v) r L K r J 

cos^ [ 1 4 , - 



I' 1 ' Si " r fi , 4 or.- ■>....,-,. ,,,-.,.2 



2tt(1 -i/) r 
sin ip 
"tt(1 - i/) ~ 



1 + -5-^ - 2^2 («x) - 2KrK 1 <Kr) - 6CV 2 )■ , (3.68) 
jl - KrK^nr) - 4CV 2 }, (3.69) 



with 



The corresponding trace is given by 

The elastic strain of this edge dislocation in a cylinder is 

E rr = r — { {1 - 2v) - + 2K 2 ( K r) + 2vktK x {kt) - (2 - 8v)Cr 2 \ 

47r(l — i/J r L K ' J 

(3.72) 

^ - 4^) ^ I 1 - ^ + ^ - ^ > ( 3 - 73 ) 



(6-8v)Cr 2 \. (3.74) 



The corresponding dilatation field reads 

7 k _ b{l-2v) siny 
2?r(l -v) r 



|l-Kr^i(Kr)-4Cr 2 |. (3.75) 



From the condition (3.38) we find for the rotation of an edge dislocation in a 
cylinder 

u z = -u = ~ i _ nrK^nr) + 4Cr 2 \, (3.76) 

ztt r I J 
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which gives an additional constant contribution to the bend-twist tensor k x in 
Eq. (3.64). The dislocation density of the edge dislocation in a cylinder is deter- 
mined by the condition (3.37) and it has the same form as the dislocation in an 
infinitely extended body (3.58). 

For convenience we give the displacement components in Cartesian coordinates 



= ^ ^(l - ktK\ (/tr)) + ^ sign (y)nrK 1 (nr) 



b f , > x 2 -y 2 ( 4 



(1 - 2i/)(lnr + K a (nr)) + *^J? (l - -±_ + 2 ^ 2 ( K r)) 



4vr(l - i/) 

+ c(x 2 (5-4^)-y 2 (l-4^))J. (3.78) 

Note that this boundary-value problem of the edge dislocation may be applied to 
describe dislocation behaviour in nanoparticles and nanowires. 



4 The relation to the nonlocal theory of dislocations 

Let me now give the relation to Eringcn's so-called nonlocal theory of dislocations [28, 
30,31]. It includes the effect of long range interatomic forces so that it can be used 
as a continuum model of the atomic lattice dynamics. Using Green's function of the 
Helmholtz equation (3.59) we may solve the field equation for every component of the 
stress field (2.46) by the help of the convolution integral: 

a ij (r)= [ G(r-r')a t3 (r')dv(r') (4.1) 
Jv 

with 

G(r -r') = \ a xz (r - r') = — K ( K (r - r')). (4.2) 

In this way we deduce Eringen's so-called nonlocal constitutive relation for a linear 
homogeneous, isotropic solid. Therefore, the gauge theoretical description of disloca- 
tions in this paper is close to the nonlocal elasticity of dislocations [28,30,31], if one 
uses the special two-dimensional Green's function (4.2). The form of the nonlocal 
kernel (dislocation density) is fixed by the incompatible distortion. Finally, we may 
conclude that the gauge theory of dislocation contains a weak nonlocality due to the 
moment stresses in the dislocation core region. This observation is in agreement with 
Kroner's opinion that the introduction of moment stresses represents the first step in 
a transition from local to a nonlocal elasticity theory [77-80] . 

The characteristic internal length in these theories is This length may be 

selected to be proportional to the lattice parameter a for a single crystal, i.e. 

KT 1 = e a, (4.3) 

where eo is a non-dimensional constant or material function which can be determined 
by one experiment [30] . For eo = we have classical elasticity and there are no effects 
of moment stresses. In this limit the nonlocal kernel must revert to the Dirac delta 
measure. In [28, 30] the choice e = 0.399 and in [32-34] the choice e = 0.25 are 
proposed (see also [2,3,38]). 
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Consequently, we can conclude that the gauge theoretical solution of the stress 
field for an edge dislocation given in this paper is also a solution in nonlocal elasticity. 
However, a difference to our approach is that Eringen's nonlocal theory works with 
compatible strain instead of incompatible strain and the strains and displacements 
have the classical form with singularities at the dislocation line. 

5 Discussion and Conclusion 

The translational gauge theory of dislocations has been employed to consider the 
straight edge dislocation in an infinitely extended body and in a cylinder. This trans- 
lational gauge theory unifies elasticity and plasticity to elastoplasticity in field the- 
oretical way. In this framework, elastoplasticity is a theory with force stresses and 
moment stresses. The size of this moment stress can be calculated even from lin- 
earized elastoplasticity theory. Using the stress function method, exact analytical 
solutions for the displacements, strain and stress fields have been derived which have 
not any artificial divergency at the dislocation core and the corresponding far fields 
agree with the classical ones for the edge dislocation. The strain and stress fields of 
an edge dislocation in an infinitely extended body are in exact agreement with those 
obtained by Gutkin and Aifantis [36, 37] through an analysis of the edge dislocation 
in the theory of gradient elasticity. Nevertheless, the gauge theoretical and the gra- 
dient theoretical approaches are different and they give slightly different results, e.g. 
the displacement field u x due to the incompatible distortion in the gauge theory of 
dislocations. On the other hand, the strain gradient elasticity used a purely compat- 
ible distortion. However, the difference cannot be too big and the numerical values 
should be of the same order. The dislocation density of an edge dislocation obtained 
in the gauge theoretical framework is in agreement with the nonlocal kernel given by 
Ari and Eringen [28-31]. From that point of view the nonlocal elasticity theory ap- 
proach is close to ours. But, on the other hand, Eringen's nonlocal elasticity contains 
compatible strains only and the strain and displacement fields have the same form as 
in a classical elasticity. In the gauge theoretical analysis the dislocation core arises 
naturally. We have calculated the localized moment stress which is the response to 
dislocation density and the far-reaching rotation gradient (bend- twist tensor) . 

The dislocation model used in this paper contains, in general, 10 material con- 
stants. The total Lagrangian of this model is given in terms of the distortion, the 
dislocation density (torsion tensor) and deWit's bend-twist tensor. By the help of 
some special choices (2.24), (2.40) and (2.43) we reduced the number of constants to 
3 material constants, namely 2 elastic constants and 1 internal length. As a result we 
obtain a dislocation theory with symmetric force stresses. The special choice (2.24) 
for the three coefficients di, a 2 and 03 in the dislocation gauge Lagrangian gives the 
core energy of the edge dislocation in an infinite medium which is in agreement with 
the expression obtained in the Peierls-Nabarro model. It is important to note that one 
can also use the total Lagrangian of this dislocation model for the gauge theoretical 
description of screw dislocations. In this turn, one reproduces all quantities of the 
screw dislocation calculated in [2]. Therefore, by using the Lagrangian proposed in 
this paper, one is able to obtain gauge theoretical solutions for the edge and screw 
dislocation which have no singularities and agree with the classical solutions in the 
far field. The (linear) gauge equation formulated in terms of the stress field is an 
inhomogeneous Helmholtz equation for the screw and the edge dislocation. 

It will be interesting to apply the nonsingular solutions of the edge and screw 
dislocation to dislocation modelling. One can use the gauge theoretical displacement 
fields to simulate the core region for HRTEM micrographs. Another interesting ap- 
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plication could be the x-ray diffraction profile analysis. One should use the materials 
tungsten, diamond or aluminium because they are nearly isotropic. In this way the 
coefficient kT 1 should be experimentally verified. Additionally it would be of interest 
to compare experimental results with the values from different theoretical approaches 
(gauge theory, strain gradient elasticity and/or nonlocal elasticity). Another possi- 
bility to determine the coefficient kT 1 should be a lattice-theoretical approach. The 
coefficient kT x expresses the interactions and the properties in the core which are 
different from those in the undeformed crystal. Consequently, the elastic moduli say 
nothing about the interactions in the core of dislocations. It should be possible to 
calculate n~ x from the coupling parameters of the lattice at the dislocation core. 
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